Is there further evidence for spatial variation of fundamental constants? 
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The detection of a spatial variation of the fine-structure constant, a, based on study of quasar 
absorption systems has recently been reported The physics that causes this a- variation should 
have other observable manifestations, and this motivates us to look for complementary astrophysical 
effects. In this paper we propose a method to test whether spatial variation of fundamental constants 
existed during the epoch of big bang nucleosynthesis. Using existing measurements of primordial 
deuterium abundance we find very weak indications that such a signature might exist, but the 
paucity of measurements precludes any firm conclusion. We also examine existing quasar absorption 
spectra data that are sensitive to variation of the electron-to-proton mass ratio fi and x = a 2 /j,g p 
for spatial variation. 



I. INTRODUCTION 

The results of a very large study using data from both 
the Keck telescope and the Very Large Telescope provide 
evidence that there is a spatial gradient in the variation of 
the fine structure constant, a = e 2 /Tic [1]. In one direc- 
tion a appears to have been smaller in the past, while 
in the other direction it appears to have been larger. 
Briefly, the method compares atomic spectra taken on 
Earth with spectra seen in absorption systems at high 
rcdshift. Any change in a results in well- understood dif- 
ferences between the two spectra. The use of many atoms 
and ions in a large number of systems ensures good con- 
trol of systematics @, Q . 

Astrophysical observations previously suggested that a 
may have been smaller in the past [4]-|7j]. However, these 
studies all used spectra taken at the Keck telescope in 
Hawaii, at a latitude of 20° N. The recently observed 
gradient in a, which we will refer to as the "Australian 
dipole" , has a declination of around —60°, which explains 
why the Keck data, restricted mainly to the northern sky, 
originally suggested a time- varying a that was smaller in 
the past. Similar studies, using a much smaller sample 
from the VLT in Chile (latitude 25° S) at first showed a 
stringent null constraint @ . More careful analysis of the 
this sample suggested that the errors should be enlarged 
by a factor of six [ljj. The recent study reported 
in [l[ is the first large-scale analysis of VLT data for co- 
variation. 

The detection of this spatial dipole motivates us to 
reinterpret, in terms of spatial variation, the existing 
studies of variation of different fundamental constants 
in other systems. In this paper we test whether fur- 
ther evidence for a dipole can be found in existing mea- 
surements of variation in electron-to-proton mass ratio, 
/i = m e /m p , and combinations of fundamental constants 
such as x = a 2 ^,g p , where g p is the proton g-factor. Fi- 
nally, we consider possible variations in the relative pri- 
mordial abundances of elements created in the during big 
bang nucleosynthesis. We develop a model for the detec- 
tion of a dipole in primordial abundances and its inter- 
pretation in terms of spatial variations of fundamental 
constants during the first few minutes after the big bang. 



II. MODEL 

We model the spatial variation of a fundamental con- 
stant X by 
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where SX/Xq = (A(r) — Ao)/Ao is the relative variation 
of A at a particular place r in the Universe (relative 
to Earth at r = 0). dx is the strength of the spatial 
variation, and E(r) describes its geometry. Note that we 
have assumed that 8X/Xq = at zero redshift, however 
this assumption should be tested for each system one 
measures using the same methods that are used at high 
redshift (e.g. by using absorbers within our own galaxy). 

Since we are fitting a dipole, we use S(r) = rcosip 
where r = ct is light-travel distance measured in giga- 
lightyears and ip is the angle between the direction of 
the measurement and the axis of the dipole. ip requires 
two parameters to specify the direction of the dipole: 
right ascension, <f>, and declination, 9, in the equatorial 
coordinates. Then cos ip is given in terms of the direction 
of the measurement and dipole axis 

H(r)=rcoflV(&,0«O ( 2 ) 
cos ip — cos(0 — 4>d) cos cos 0d + sin 9 sin 9d ■ 

We have additionally assumed that the effect of any spa- 
tial variation of fundamental constants increases with 
distance. This is model dependent for large redshifts: 
we use the standard Acdm cosmology parametrized by 
WMAP5 [ll| to determine comoving distance r. In these 
coordinates the Australian dipole of a- variation found by 

Qis 

Set 

— = (1.10 ±0.25) x l(T 6 r cos ip , (3) 
a 

where ip is the angle between the direction of the 
measurement and the axis of the Australian dipole, 
(17.4 (0.6) h, —58 (6)°) in equatorial coordinates. 

For any given set of data measuring variation of funda- 
mental constants, one should ask whether a dipole model 
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fits the data better than a monopole or a null hypothe- 
sis (that is, no variation). However, our goal here is to 
see if any of the existing data provide further evidence 
for the spatial variation of fundamental constants seen 
in [l[ . Therefore it is also reasonable to question whether 
a dipole model with the axis specified by the Australian 
dipole provides a better fit to the data. There are some 
good theoretical justifications for such a procedure. For 
instance, the constants may vary because they are cou- 
pled to a (dimcnsionless) scalar field $ which varies over 
space-time — for example, the quintessence field <£>/c 2 or 
a dimcnsionless dilaton field. In this case the axis of the 
dipole is the direction of its gradient V<&, and a funda- 
mental constant X is coupled to its variation via 

S -^ = k x 6<P, (4) 

where kx is a dimensionless coupling coefficient. Our 
dipole model now requires (S$(r) ~ S(r) but all constants 
will vary in the same direction (i.e. along the Australian 
dipole) . 

In this paper we discuss the constants a, [i, and the di- 
mensionless mass ratio X q = m q / Aqcd, where m q is the 
light-current quark mass and Aqcd is the position of the 
Landau pole in the logarithm of the running strong cou- 
pling constant, a s ~ 1/ ln(AQCD'Vic). In the Standard 
Model the electron and quark masses are proportional to 
the vacuum expectation of the Higgs field, while the pro- 
ton mass m p is proportional to Aqcd (if we neglect the 
~ 10% contribution of the quark masses). Relative vari- 
ation of X q is then approximately equal to the relative 
variation of m q /m p and fi — m e /m p . We can relate the 
relative variation of different constants by equations like 

k^ = R k a , (5) 

where the Rx can be determined from observations and 
compared with theories of spatial variation. 



III. H 2 ABSORPTION LINES 

Limits on variation of /i at high redshift come from 
comparison of different rotational-electronic excitations 
in molecular hydrogen (H2). The most recent determi- 
nations use the laboratory wavelengths and sensitivity 
coefficients presented in |L^ - [l4| to determine 5/J.//J, in four 
different quasar absorption systems fl5l . [l6j . 

To determine which model of variation provides the 
best fit to the Sfx/fj, data, we can use x 2 l v as a goodness- 
of-fit parameter. Here v is the number of degrees of free- 
dom; for four data points (quasar absorption systems), 
v = 4 — p where p is the number of parameters in the 
model. For example, fitting a monopole to the data re- 
duces x 2 l v from 1.21 (null hypothesis) to 1.08 (one pa- 
rameter), which explains why the data mildly supports a 
non-zero "detection" oiSfi/fj, = 3.4(2.7) x 10~ 6 . 



Alternatively, we can test whether the data supports 
a dipole model in the direction of the Australian dipole. 
Using the model (pj with the dipole direction fixed at 
(17.4 h, -58°) we obtain d M = 2.6 (1.3) x 10 _6 /Gyr, with 
a reduced x 2 /v = 0.33. We can compare this to the null 
hypothesis using an F-test, which allows us to find the 
probability that x 2 has improved by chance. We define 
Fhy 

F(p 2 > Pi = 277 \ > 6 

where p\ and pi are the number of parameters in the re- 
spective models and n is the number of data points (in 
the present case, pi = 0, P2 = 1, and n = 4). Note that 
we cannot use this test to determine whether the dipole 
fit is better than the monopole (although x 2 l v suggests 
that it is), since the models compared using the F-test 
must be "nested" : it must be possible to generate model 
1 from model 2. In any case, since there is no detec- 
tion the null hypothesis is really the better comparison. 
We obtain F = 11.6 and the probability that this result 
occurs due to chance is 4%. 

Clearly the Australian dipole provides a good fit to 
the data, but what does the data alone say? To find 
the direction preferred by the data, we use the three 
parameter fit (Q]); the result is rf M = 3.3(1.5) x 10~ 6 , 
4> d = 16.7(1.5) h, 9 d = -62(5)°. This direction is con- 
sistent with the Australian dipole. With only one de- 
gree of freedom remaining, x 2 = 0.02. The errors for 
each model parameter are found by varying the param- 
eter locally near the best-fit value until x 2 = Xmin + 1j 
while keeping all other parameters fixed at their respec- 
tive best-fit values. This ignores correlations between 
the parameters, which must be important; for example if 
dp = (two standard deviations from its best-fit value) 
then cj)d and 6d are completely unspecified. 

While the preceding paragraphs could be interpreted 
as evidence for spatial variation of /x, a few things must 
be noted. Most obviously, we have a very small sample 
of data, which in principle can be improved although 
there are significant observational challenges in obtaining 
data of sufficient quality at high redshift [ill E3 • Also 
of importance to this work is the lack of sky-coverage 
represented in the sample. All four sources are below 
the equator, and three of them are within 25° of each 
other. Therefore many more data points over a wider 
sky coverage are needed to obtain a significant detection. 

Once a detection is confirmed, we will be able to ex- 
tract R^] of Eq. ([5]). For example, with the suggested 
by this work and d a given by ©, R* = 2.6 x lO" 6 / 1 - 1 x 
10" 6 = 2.4. However, in Section [V] we discuss stringent 
limits from studies of ammonia inversion spectra which, if 
confirmed, place very strong limits on cosmological vari- 
ation of p and hence R^ . 
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IV. 21-CM RADIO VS. UV LINES 

The variation of a combination of fundamental con- 
stants is probed by comparing neutral hydrogen (HI) 21- 
cm radio and UV metal lines that have been redshiftcd 
into the optical band (see 0, Qj|). The ratio of these 
transition frequencies is proportional to x = o? [ig v . Only 
a handful of quasars have been found for which both op- 
tical and radio absorption are observed; the nine systems 
with good data are presented in [l9j]. These span the 
redshift range z = 0.23 to 2.35. 

It was found in [l9| that there is far more scatter in the 
data than would be expected from the statistical errors: 
the best monopole (single parameter) fit gives x 2 l v = 8-1 
(x 2 = 64.6), where ~ 1 would be expected if the er- 
rors were Gaussian and the model correct. Therefore [l9| 
concluded that there are large systematic errors, which 
they trace to the assumption that the 21-cm absorption 
component is at the same physical location as the UV 
component. It is supposed that the scatter seen in the 
data mainly comes from this assumption being incorrect; 
that is, the 21-cm and UV absorbing gases can be ran- 
domly offset in space and have quite different velocities 
(although see [20l, l2l|j). 

Having enlarged the errors in order to force x 2 l v = 1 
for the monopole model, it is quite unsurprising that [l9[ 
found no evidence for a dipole: the monopole model has 
been made to fit as well as could possibly be expected. 
Therefore we have reanalysed the data without the addi- 
tional systematic to test whether any significant dipole 
exists. Firstly, using the axis of the Australian dipole, we 
find x 2 reduces from 110 (zero parameters) to 101 (one 
parameter), x 2 l v is actually worse than for the monopole 
fit. Fitting the three parameter dipole ([T]) provides part 
of the explanation for why the Australian dipole is such 
a poor fit: the best fitting dipole is at (7.3 h, 23°), an 
angle of 140° from the Australian dipole. For this model 
X 2 lv — 5.5, which shows that even with the best-fit 
dipole systematics still dominate. 

We conclude that the detection of a dipole in x of real 
significance will require much more data than the other 
systems discussed in this paper, due to the need to ac- 
count for the large systematic errors caused by the veloc- 
ity offset between the HI and metal absorption lines. 



V. INTERPRETING SINGLE MEASUREMENTS 

In this section we briefly discuss how individual mea- 
surements of combinations of fundamental constants can 
be compared with the dipole model presented in Sec- 
tion [TTJ When there are only one or two samples, it is 
not possible to determine whether a dipole fit to the data 
is appropriate. However, one can still compare the mea- 
surement with the variation that the Australian dipole 
would lead us to expect. 

Two new measurements comparing HI 21-cm with 
neutral carbon (CI) absorption lines at z = 1.36 and 



z = 1.56, along lines of sight to quasars Q2237-011 
and Q0458-020, have been reported recently [22[ . As 
discussed in Section [TV] these measurements probe the 
combination of constants a 2 [ig p . In Table UJ we compare 
the results of these measurements with the variation ex- 
pected according to the Australian dipole and the model 
presented in Section [TTJ To obtain the model prediction, 
we calculate the expected variation of a according to ([3]) 
and leave the factors Rg p and R^, defined by equations 
such as (O, for later determination. The type of anal- 
ysis presented here could eventually yield values for the 
ratios Rg p and i?^, but a larger sample of measurements 
is required. 

Measurements of /z-variation have been made by com- 
paring inversion lines of ammonia, NH3, with rotational 
lines of CO, HCO + , and HCN molecules. The inver- 
sion spectrum of ammonia has an enhanced sensitivity 
to variation of /i because it depends on the exponen- 
tially small probability of tunneling of the three hydrogen 
atoms through the potential barrier [27J • In Table UJ we 
present two measurements of /x variation based on this 
method. 

Two other measurements are shown in Table [JJ The 
combination (a 2 j /j,) 157 g p can be extracted by compar- 
ison of OH 18-cm and HI 21-cm lines; this has been 
performed for the grav itational lens toward PMN J0134- 
0931 at z = 0.765 |25j| . A measurement of (a 2 /M) 185 3 P , 



derived from comparison of conjugate-satellite OH 18- 
cm lines, has yielded a 2.6<r non-zero detection towards 
PKS 1413+135 26]. We note that this absorber lies at 
~ 81° to the Australian dipole and is at low redshift 
(r ~ 2.9 Gyr), therefore minimal variation in fundamen- 
tal constants is expected. 



VI. BIG BANG NUCLEOSYNTHESIS 

In principle one can search for spatial variation in the 
primordial abundance of any element produced during 
BBN (the most common are D, 3 H, 3 He, 4 He, 7 Li). The 
best limits on primordial deuterium abundance are de- 
rived from isotope-shifted Ly-a spectra in quasar ab- 
sorption systems, where the absorbers are at z ~ 2.5. 
By contrast 4 He is observed in ionized (HII) regions of 
low-metallicity dwarf galaxies at z < 0.01, while the pri- 
mordial 7 Li abundance is determined from metal-poor 
Population II stars in our galaxy. Therefore, in practice 
only the deuterium abundance is currently measured at 
large enough rcdshifts to expect to see any effect of spa- 
tial variation of fundamental constants. 

In the case of primordial abundances we need to ex- 
tract an average value from the data (monopole term) as 
well as any potential dipole. Therefore we use the model 

log 10 (a/H) = m a + d a 3(r) (7) 

= m a + d a r cos ip{<f>d,6d) , 

where log 10 (a/H) is the primordial abundance of element 
a relative to hydrogen abundance at a particular place 
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TABLE I. Comparison of expected variation, given by Eq. ([3]), and measured variation of fundamental constants in different 
systems. Each measurement corresponds to a single absorption system. Rg p and are defined by equations like ([5|. 



System 


Constant 


Expected variation 
(xlO" 6 ) 


Measurement 
(xlO" 6 ) 


Ref. 


HI 21-cm + CI 

NH3 inversion 

HI 21-cm + OH 18-cm 
OH 18-cm 


/'• 


1.12 x {2 + R° + R£ r ) 
-5.20 x (2 + R% + R« p ) 
-5.47 B° 
1.34 R« 

-1.04 x (3.14 - 1.57/^+7?^) 
0.50 x (3.70 - 1.85/?° + R* p ) 


6.64 ± 0.84 stat ± 6.7 sys 
7.0 ± I.8 Btat ± 6.7 sys 

< 1.8 x 10" 6 (2a) 

< 1.4 x 10~ 6 (3cr) 
4.4 ± 3.6 stat ± 10 sys 

-11.8 ±4.6 


u 

[23] 
[24] 
[25] 
[26] 



r, m a is the average primordial abundance, d a is the 
strength of the spatial variation, and 3(r) describes the 
geometry of the spatial variation @. 

In the future one might hope to measure primordial 
abundances of the other elements at high redshift. The 
production of each of these elements has a different sensi- 
tivity to fundamental constants, therefore a complete set 
of data can, in principle, simultaneously measure sev- 
eral fundamental constants at the time of nucleosynthe- 
sis. The dependence of primordial abundances on funda- 
mental constants such as a and m q / Aqcb (the ratio of 
light-quark mass to the pole in the running strong cou- 
pling constant) is the subject of current research [2Sl - 43lj | . 
An observed spatial variation in primordial abundance 
(non-zero d a ) can be related to the variation of a funda- 
mental constant X at that position in space at the time 
of big bang nucleosynthesis using the relationships 



d a = 



flbg 10 (o/H) 91og 10 (a/H) d\nX 



9 In A 



dZ 



(8) 



where <91og 10 (a/H)/cnn X is determined from theory. 

We have performed a preliminary search for such a 
dipole using the seven best determinations of primor- 
dial deuterium abundance presented in [32[. It is known 
that the standard (monopole) model gives a minimum 
X 2 of 19.13 rather than the ~ 6 that would be expected 
if the errors were purely statistical, therefore there is at 
least a possibility that the dipole model of spatial vari- 
ation could produce a better fit to the data. As in the 
cases presented in earlier sections, we can investigate this 
possibility using \ 2 jv as a measure of the fitness of our 
model. The results are shown in Table HU We find that 
the dipole model with direction fixed by the Australian 
dipole [U is not significantly preferred over the simple 
monopole model, with x 2 l v increasing very slightly. On 
the other hand, if we let the data choose the dipole direc- 
tion, in effect solving for m, d, 4>d and 9d simultaneously, 
we only have three remaining degrees of freedom, and 
X 2 lv increases substantially. Interestingly, the dipole is 
found to point in the same direction (within errors) as 
the Australian dipole (see Table ITT]) . 

In Fig. [T] we show the measured abundances against 
rcosifj- The dipole fit (solid line) fits better than the 
monopole fit (dashed line) . We clearly need more sky cov- 
erage: all of the data points lie at angles 55° < -0 < 125° 



FIG. 1. Deuterium abundance vs. r cos ip where ip is the 
angle between the measurement direction and the dipole of 
Ref. [J| . Dashed line: monopole fit; solid line: dipole fit (fixed 
direction). Fit parameters are given in Table iHl 
-4.3 
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from the dipole axis. The distances are similar for all the 
data points (r ~ 11 Gyr), therefore our results do not 
vary significantly if we remove the distance dependence 
and let S(r) = cos-0- We note in passing that if the two 
outlying points in Fig. Q] are removed and only five re- 
main, x 2 l v — 2-00 for the monopole fit and 0.44 for the 
dipole fit: m = -4.55 (2) and d = 0.010 (4). 

A more reasonable approach to the x 2 discrepancy is 
to increase the error bars to account for an unknown 
systematic. We have done this by adding a constant 
error er sys to all the statistical errors in quadrature: 
a 2 ot = <Js tat + o~ 2 ys . Increasing the errors such that the 
dipole fit has x 2 = 5 requires an additional systematic of 
cr sys = 0.11 to be added to all the data points. An in- 
crease in the uncertainty of this magnitude destroys any 
justification for fitting a dipole to the data. 

Finally, we interpret our "detection" of a dipole in the 
primordial deuterium abundance, aligned with the Aus- 
tralian dipole, in terms of variations in fundamental con- 
stants. The parameter dp = 0.0045 (35) and we may 
interpret this in terms of variation of a. Using the result 
<91n(D/H)/<91na = 3.6 30] in Eq. © we obtain 



da 



0.0045 



a 3.6/ In 10 



3(r) = 0.003 3(r) . 
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TABLE II. Calculated best fits to the primordial deuterium abundance data using various models. The model parameters, 
including mn, dr> and the direction of the dipole (right ascension and declination) from (flj, as well as minimum y 2 and \ 2 
per degree of freedom (y) are shown. The monopole model is the standard interpretation (mo only). The dipole model was 
tested using fixed and varying dipole directions. The fixed dipole direction corresponds to the best fit of [ij: (17.4 h, —58°) in 
equatorial coordinates. Errors in each parameter are found ignoring correlations; see Section iHll text for details. 





rno 


(Id 


R.A. 

(hours) 


Decl. 
(deg) 


x 2 


x> 


Monopole 


-4.55 (2) 








19.13 


3.19 


Dipole (fixed direction) 


-4.55 (2) 


0.0045 (35) 






17.50 


3.50 


Dipole 


-4.56(2) 


0.0054 (29) 


15.5(1.6) h 


-14(51)° 


15.73 


5.24 



In other words, a system 1 Gyr in the direction of the 
dipole axis had a fine-structure constant different from 
the present-day laboratory value by da/a = 0.003 at the 
time of big bang nucleosynthesis. 

This effect is much larger than that seen directly in 
quasar absorption spectra: 5a/ a ~ 10~ 6 S in Yet 
even if this signal is confirmed, there is no inconsistency 
here. The direct measurements of [l[ probe physics oc- 
curring a few billion years ago, while the primordial abun- 
dance measurements probe the first few minutes after the 
big bang. 

We chose the fine-structure constant since spatial vari- 
ation of a has been reported. However big bang nu- 
cleosynthesis is more sensitive to X q = to 9 /Aqcd, the 
dimensionless ratio of light-quark mass to the pole in 
the running coupling constant, and it is generally more 
reasonable to interpret the variation in terms of smaller 
shifts in these constants. For example using the result 
<91n(D/H)/<91nJTg = 7.7 we obtain 

^2 = 0.0013 (10) 5(r). 

Clearly the detection and confirmation of spatial vari- 
ation in fundamental constants at the time of big bang 
nucleosynthesis requires a large improvement in the num- 
ber of deuterium measurements and their coverage of the 
sky. However, it may also be possible to detect other el- 
ements produced in BBN such as 3 He, 4 He, 6 Li, and 7 Li 
in high-redshift quasar absorption spectra and determine 
their primordial abundances. If the response of BBN to 
variation of constants is known, one can readily predict 
the spatial variation of their abundances given a model 
of fundamental constant variation. We reproduce some 
of the known response functions from the theory [29|, l30j 
in Table EH 

We can see, for example, that a variation of 1% in 
a would amount to a 1.9% change in Y p (the mass ra- 
tio 4 He/H - for all other elements the number ratio is 
presented) and a —11% change in 7 Li abundance. This 
provides an independent means of verifying a particu- 
lar change in fundamental constants. The primordial 
7 Li abundance evidently provides a particularly sensitive 
probe of fundamental constants. In Table Mil we also 
present the sensitivity of primordial abundances to the 
baryon-to-photon ratio r\ (relative to the WMAP value 



TABLE III. The sensitivity of relative variation in primor- 
dial abundances to relative variation in fundamental constants 
dlnK/dlnX. Here Y a = (D/H, 3 He/H, Y p , 6 Li/H, 7 Li/H) 
are number ratios of primordial isotope abundances to hydro- 
gen, except Y p which is the mass ratio 4 He/H. 



X 


D/H 


3 He/H 


Y P 


6 Li/H 


7 Li/H 




a 


3.6 


0.95 


1.9 


6.6 


-11 


[30] 


x q 


7.7 




-0.95 




-50 


[29] 


V 


-1.6 


-0.57 


0.04 


-1.5 


2.1 


[30] 


G 


0.94 


0.33 


0.36 


1.4 


-0.72 


[30] 



?7 = 6.1 x 10 10 ), which prior to WMAP was constrained 
chiefly by BBN data. 



VII. CONCLUSION 

The existing H2 data from quasar absorption spectra 
shows hints that there may be a dipole in /i-variation 
with an axis corresponding to that of the gradient in 
a found in [l|. On the other hand x — a 2 [ig v varia- 
tion, taken from HI 21-cm data, has a best-fit dipole 
whose axis does not correspond to that of the Australian 
dipole, although in this case systematics heavily domi- 
nate. We have demonstrated that it is possible to in- 
fer spatial variation of fundamental constants during big 
bang nucleosynthesis from high-redshift measurements of 
primordial abundances. Although the existing deuterium 
data does not necessarily support the dipole interpreta- 
tion, with the significance of the dipole model being sim- 
ilar to that of the monopole model, it is interesting that 
the preferred axis is consistent with the direction of the 
Australian dipole. There is a strong impetus now to per- 
form measurements of relative primordial abundance at 
high redshifts of as many elements as possible in as many 
different spatial directions as possible. 

Finally, we note that it may be possible to observe a 
spatial dipole in other cosmological systems [33[ . For ex- 
ample, a-variation may be seen in the CMB anisotropy 
if a high-enough sensitivity can be reached. Although 
the results of [l| (interpreted as strictly spatial variation) 
suggest that accuracy at the level 10 -6 will be required, 
if the hints from BBN turn out to be real then there is 
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an additional redshift (time)-dependence that could in- 
crease the variation at the time of the CMB substantially. 
Another possibility is that if the observed a-variation is 
related to the cosmological constant, and hence the ac- 
celerated expansion of the Universe, it may be possible 



to see a dipole in the redshift-luminosity relationships of 
Snla supernovae data (see, e.g. [34j|). 

We thank M. T. Murphy for useful discussions. This 
work is partly supported by the Australian Research 
Council. 
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